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Abstract The future Compact Linear Collider (CLIC) of¬ 
fers a possibility for a rich precision physics programme, in 
particular in the Higgs sector through the energy staging. This 
is the first paper addressing the measurement of the Standard 
Model Higgs boson decay into two muons at 1.4 TeV CLIC. 
With respect to similar studies at future linear colliders, this 
paper includes several novel contributions to the statistical 
uncertainty of the measurement. The later includes the Equiv¬ 
alent Photon Approximation and realistic forward electron 
tagging based on energy deposition maps in the forward 
calorimeters, as well as several processes with the Beam- 
strahlung photons that results in irreducible contribution to 
the signal. In addition, coincidence of the Bhabha scattering 
with the signal and background processes is considered, al¬ 
tering the signal selection efficiency. The study is performed 
using a fully simulated CLIC_ILD detector model. It is shown 
that the branching ratio for the Higgs decay into a pair of 
muons BR(// —> /i jl ) times the Higgs production cross- 
section in IT IT-fu si oil o(Hvv) can be measured with 38% 
statistical accuracy at \fs =1.4 TeV, assuming an integrated 
luminosity of 1.5 ab 1 with unpolarised beams. If 80% elec¬ 
tron beam polarisation is considered, the statistical uncer¬ 
tainty of the measurement is reduced to 25%. Systematic 
uncertainties are negligible in comparison to the statistical 
uncertainty. 

1 Introduction 

Measurements of Higgs branching ratios, and consequently 
Higgs couplings, provide a strong test of the Standard Model 
(SM) and possible physics beyond. Models that could possi¬ 
bly extend the SM Higgs sector (Two Higgs Doublet model, 

*This work was carried out in the framework of the CLICdp col¬ 
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Little Higgs models or Compositeness models) will require 
Higgs couplings to electroweak bosons and Higgs-fermion 
Yukawa couplings (coupling-mass linearity) to deviate from 
the SM predictions [I]ED. 

The Compact Linear Collider (CLIC) represents an ex¬ 
cellent environment to study properties of the Higgs bo¬ 
son, including its couplings, with a very high precision (3] 
a. Measurements of rare H —> fi 1 jl decays are particu¬ 
larly challenging because of the very low branching ratio 
of 2xl0“ 4 predicted in the SM Q for a Higgs mass of 
126 GeV. Current results indicate that the LHC was not able 
to access Higgs coupling to muons (gH/j^i), based on the runs 
at 7 TeV and 8 TeV centre-of-mass (CM) energies |6|. Pro¬ 
jections for the HL-LHC, assuming 300 fb 1 and 3 ab 1 of 
data, predict uncertainties of 23% and 8% respectively for 
the gHuii coupling (7). In order to provide the best physics 
reach in the shortest time and for an optimal cost, the opera¬ 
tion of the CLIC accelerator is foreseen in energy stages of 
350 GeV, 1.4 TeV and 3 TeV ©. At 1.4 TeV and 3 TeV, suffi¬ 
ciently large Higgs boson samples can be produced to allow 
studies of rare Higgs decays . A sample of 3.7 x 10 5 Higgs 
bosons can be produced at 1.4 TeV CM energy, for an inte¬ 
grated luminosity of 1.5 ab -1 with unpolarised beams. With 
the expected instantaneous luminosity of 3.2x 10 34 cm 2 s 1 
this can be achieved in approximately five years of detector 
operation, with 200 running days per year and an effective 
up-time of 50% 0. The signal sample size will be doubled 
at 3 TeV CM energy due to rising cross-section for WW- 
fusion flOl . 

A similar study has been performed at 3 TeV CM energy 
(TO) . Compared to the study at 3 TeV, several challenges for 
the measurement of H —> j-L'j-l at CLIC are discussed for 
the first time in this paper. First, background processes with 
photons in the initial state simulated using both the expected 
Beamstrahlung spectrum at CLIC and the Equivalent Photon 
Approximation (EPA) CUE], were considered. Forward 
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electron tagging (Section[5]i leads to a rejection of 48% and 
42% of the e~e + —» e~e + p + p “ and e ± y —> e ± p + p^ back¬ 
ground events, respectively. The impact of Bhabha scatter¬ 
ing events on the rejection of events with forward electrons 
is investigated. 

The paper is organised as follows. In Section [2]the sim¬ 
ulation tools used for the analysis are listed and in Section 
[3] the CLIC_ILD detector model is briefly described. Signal 
and background processes and event samples are discussed 
in Section [4] Tagging of background high-energy electrons 
is described in Section [5] Event preselection and the final 
selection based on a Multivariate Analysis (MVA) approach 
are described in Section [6] The di-muon invariant mass fit 
and the extraction of the statistical uncertainty of the mea¬ 
surement are described in Section [7] In Section 0 the im¬ 
pact of electron polarisation on the statistical uncertainty of 
the o(Hvv) x BR (H —> pp ) measurement is described. 
Systematic uncertainties are discussed in Section[9] followed 


by the conclusions in Section 10 


2 Simulation and analysis tools 

Higgs production through WW -fusion is simulated in Whizard 
1.95 mm including a realistic CLIC beam spectrum and 
initial state radiation. The generator Pythia 6.4 ca is used 
to simulate the Higgs decay into two muons. The CLIC lu¬ 
minosity spectrum and beam-induced processes are obtained 
by GuineaPig 1.4.4 da. Background events are also gener¬ 
ated with Whizard using Pythia 6.4 to simulate the hadroniza- 
tion and fragmentation processes. Simulation of tau decays 
is done by Tauola fT71 . The CLIC_ILD detector simulation 
is performed using Mokka fl8l based on Geant4 Q9}. Be¬ 
fore digitisation of the detector signals, pile-up from yy —► 
hadrons interactions is overlaid on the physics events. The 
particle flow algorithm, PandoraPFA |[20l [2Tj| is employed 
in reconstruction of the final-state particles within the Mar¬ 
lin reconstruction framework ll22l . The TMVA package ll23l 
is used to separate signal from background by MVA of sig¬ 
nal and background kinematic properties. 


cover polar angles down to 7°. They have a single point res¬ 
olution of 7 pm, and together with the TPC improve the 
tracking accuracy in the r(j) plane. In order to provide preci¬ 
sion tracking and vertexing closer to the beam-pipe, a Ver¬ 
tex Detector capable of an impact parameter resolution of 
3 pm l26l is foreseen. Calorimetry at CLIC is based on 
fine-grained sandwich calorimeters optimized for particle- 
flow analysis (PFA). PFA is based on reconstruction of four- 
vectors of visible particles, combining the information from 
precise tracking with highly granular calorimetry. The de¬ 
tector comprises a central solenoid magnet, with a field of 4 T. 
High muon reconstruction efficiency of 99%, for muons above 
7.5 GeV, is achieved by combining information from the 
central tracker (TPC plus silicon tracker) with information 
provided by the iron yoke instrumented with the 9 layers of 
Resistive Plate Chamber detectors. 

In principle, hadrons produced in the interaction of the 
beam-induced photons affect the TPC occupancy and conse¬ 
quently the muon reconstruction efficiency. However, in the 
studied sample of muons from H —> p p decays, muon 
reconstruction efficiency is above 99% in the barrel region, 
in the presence of yy —> hadrons. 

The average muon transverse momentum resolution for 
the signal sample is A (1/ pj) = 3.3 x 10 _5 GeV _1 in the bar¬ 
rel region. The impact of transverse momentum resolution 
on the statistical uncertainty of a(Hvv) x BR (H —> p p ~) 
measurement is discussed in Section[9] 

In the very forward region of the CLIC_ILD detector, 
below 0=8°, no tracking information is available. The re¬ 
gion between 0.6° and 6.3° is instrumented with the two 
silicon-tungsten sampling calorimeters, LumiCal and Beam- 
Cal 1271 . for the luminosity measurement, beam-parameter 
control, as well as for the tagging of high-energy electrons 
escaping the main detector at low angles. Together with the 
very forward segments of the electromagnetic calorimeter 
covering the polar-angle region between 6.3°and 8°, it is 
possible to suppress the four-fermion SM background with 
the characteristic low-angle electron signature. The simu¬ 
lation of the very-forward electron tagging is described in 
Section^ 


3 The CLIC ILI) detector model 

The ILD detector concept l24l is modified for CLIC accord¬ 
ing to the specific experimental conditions at higher ener¬ 
gies 0. The subsystems of particular relevance for the pre¬ 
sented analysis are discussed here. A complete description 
of the CLIC_ILD detector can be found in lf25l . 

The main tracking device of CLIC_ILD is the Time Pro¬ 
jection Chamber (TPC) providing a point resolution in the 
r(j) plane better than 100 pm, for charged particles in the 
detector angular acceptance 0. Additional silicon trackers 


4 Event samples 

At sjs =1.4 TeV the SM Higgs boson is predominantly pro¬ 
duced via WW -fusion (FigurefTJ. The effective cross-section 
for Higgs production in WW-fusion is 244 fb without beam 
polarization. The Higgs production cross-section above 1 
TeV can be measured with a statistical precision better than 
1% as shown in 0. The e + e~ —> Hvv,H —> p + p~ signal 
statistics are expected to be small (of the order of a few tens 
of events) because of the small branching fraction for this 
particular decay. 
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Fig. 1 Feynman diagram of the Higgs production in VFVF-fusion and 
the subsequent Higgs boson decay to a pair of muons. 


V e 



We have simulated a sample of 24000 signal events, rou¬ 
ghly corresponding to 300 times the number of events ex¬ 
pected in 1.5 atV 1 of data. This is needed in order to pro¬ 
vide an adequate description of the signal Probability Den¬ 
sity Function (PDF) (Section [7TTj ). The signal and the domi¬ 
nant background processes are listed in Table |T] For each of 
the background processes, samples of 2 ab 1 are generated. 


In addition to the processes listed in Table |T] we have 
considered s-channel e + e~ —> ii ji production, as well as 
several processes with tau pair in the final state e+e~ ->T+T“, 
e + e~ v t v t t + t - , e + e~ —► e“e + T + T~. Tau decays be¬ 
come relevant if both taus decay into two muons which hap¬ 
pens in ~3% of cases (28l . However, the invariant mass of 
the di-muon system will not match the Higgs mass window 
considered in this analysis (see Section 6.11. The same holds 
for e + e~ —> p 1 /i production. 


In addition, processes with quasi-real photons in the ini¬ 
tial state are simulated using the EPA. In this analysis, such 
events are grouped together with the processes involving 
Beamstrahlung photons. In this way, processes with roughly 
similar kinematic characteristics are grouped together, as sh¬ 
own in Table [l] The notation e + y represents the sum of 
cross-sections for the processes with either e~y or e + y in 
the inital state. 


At y/s = 1.4 TeV, the Higgs boson is also produced via 
ZZ-fusion, with a cross-section of about 10% of the Higgs 
production cross-section in WW -fusion. However, on a test 
sample of 300 ZZ-fusion events followed by the Higgs de¬ 
cay to a pair of muons, not a single event passed the selection 
described in Sections 6.1 and 6.2 This implies an efficiency 
smaller than 1.2% (95% CL) for this channel equivalent to 
less than 0.1 events passing the final selection. Therefore, 
the Higgs production through ZZ-fusion is not considered 
relevant for this analysis. 


Photons, dominantly emitted by Beamstrahlung, produce 
incoherent pairs deposited mainly in the low-angle calorime¬ 
ters. On average, 1.3 two-photon interactions producing had¬ 
ronic final states occur per bunch crossing (29) which may 
affect the muon reconstruction in the tracking detectors. These 
hadrons are included in the analysis by overlaying 60 bunch 
crossings in the simulation, before the digitisation and event 


reconstruction phase. These events, as well as other physics 
events, are passed trough the full detector simulation |[30l . 


Table 1 List of considered processes with corresponding cross- 
sections. Cross-section values marked by * are generated with the 
additional kinematic requirements: 100 GeV< m ^ <150 GeV, and 
8°<0j 1 <172°, where stands for di-muon invariant mass and 8^ 
is the polar angle of the reconstructed muon. The cross-sections for 
all processes with photons in the initial state include both Beam¬ 
strahlung and processes with EPA photons. Cross-sections for pro¬ 
cesses e ± y —> e ± /r + jt l~ and e±y —> p + p~ represent the sum 

of cross-sections for the processes with both initial states e~ y and e + y. 


Process 

<y(fb) 

e + e~ —> Hvv.H —> p+p~ 

0.0522 

e + e~ —> vvp+p~ 

129 


24.5* 

e ± y^i 

1098* 

e ± y->e ± v fl v, i ^+p- 

30 

yy^ v^p+p- 

162 

e + e —> e + e v^v^p+p 

1.6 


5 Tagging of EM showers in the very forward region 

In the polar angle region below 0=8°, tracking information 
and hadronic calorimetry are not available. The four-fermion 
background e + e~ —> e~e jl of multiperipheral type and 

similar processes like e ± y—> can fake the miss¬ 

ing energy signature of the signal if the final state electrons 
(spectators), emitted at the polar angles smaller than 0=8°, 
escape undetected. 

Electron detection in the very forward region involves 
the reconstruction of electromagnetic showers in the pres¬ 
ence of intense beam-induced background depositing in the 
very forward calorimeters a large number of low-energy par¬ 
ticles, mostly incoherent pairs from Beamstrahlung 13T1. This 
deposition amounts to several hundred thousand of e + e~ 
pairs per bunch crossing (32). 

Furthermore, Bhabha events where one or both electrons 
are detected in the very forward calorimeters may occur in 
coincidence with either signal or background, even within 
the 10 ns time stamp. Tagging of such Bhabha electrons will 
result in the rejection of signal (or background). In order to 
prevent significant loss of signal statistics, the electron tag¬ 
ging was optimized to identify showers with energy higher 
than 200 GeV and a polar angle above 1.7° onlyQ Under 
these requirements, the loss of the number of signal events 
due to tagging of Bhabha electrons amounts to 7%. Out of 
these 7%, in slightly more than a half of events one elec¬ 
tron is added to the final state and, in the remainder two 

’Simulation with the reconstruction algorithm from Ref. ED 
shows that assuming these cuts, the reconstruction efficiency, in Beam- 
Cal and LumiCal is above 98%, with a negligible fake rate. 
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Bhabha electrons are added. Table [2] shows rejection rates 
for signal and background obtained by the forward elec¬ 
tron tagging due to Bhabha pile-up. 

In conclusion, very forward tagging of high-energy elec¬ 
trons serves to half the fraction of background with spectator 
electrons, with a moderate loss of signal of 7% in the pres¬ 
ence of Bhabha coincidence. 


Table 2 Rejection rates for signal and background by the forward 
electron tagging. 


Process 

Rejection rate 

e~e + —> e~e + p + p~ 

48% 

e ± y —> e+p+p- 

42% 

Signal 

7% 


6 Event selection 

The event selection is done in two steps. First, a preselec¬ 
tion is performed aiming to suppress background originating 
from beamstrahlung as well as the processes with spectator 
electrons described in Section [5] The final event selection 
uses a multivariate classifier based on boosted decision trees 
(BDT) to suppress remain background processes on the ba¬ 
sis of their kinematic properties. 

6.1 Preselection 

In order to suppress the impact of the beam-induced back¬ 
ground, only reconstructed particles with transverse momenta 
pj > 5 GeV are used in the analysis. Furthermore, the pre¬ 
selection of events was made by requiring a reconstruction 
of exactly two muons in the event, with an invariant mass 
of the di-muon system in the window centered around the 
Higgs mass 105-145 GeV. In addition, the absence of tagged 
electrons with energy above 200 GeV and polar angle above 
1.7° is required in order to suppress background with spec¬ 
tator electons emitted in the very forward region of the de¬ 
tector. 


6.2 MVA selection 

As a second step in the event selection, MVA techniques are 
used based on the BDT classifier implemented in the TMVA 
package. From the signal sample, quater of all events are re¬ 
served for TMVA training, as well as 0.5 ab 1 of each back¬ 
ground. The following observables were used for the clas¬ 
sification of events, similar to the CLIC study at yA = 3 TeV 

03 : 


- visible energy of the event excluding the energy of the 
di-muon system, E v j s , 

- transverse momentum of the di-muon system, pT(fifi), 

- scalar sum of the transverse momenta of the two selected 
muons, pr(ni) + Pt{H 2), 

- boost of the di-muon system, j3 M|U = \ppp \ /E^, 

- polar angle of the di-muon system, 

- cosine of the helicity angle, cos G*. 

The process e + e~ —> vvf.l 'j.1 , with the same final state 
as the signal, represents an irreducible background and can 
not be substantially suppressed before the invariant mass fit 
of the di-muon system. The process yy —> V fl V ll p + /l~ has 
a similar final state, but a different CM energy distribution 
in the initial state, since it involves Beamstrahlung or EPA 
photons rather than initial electrons. This leads to a differ¬ 
ent distribution of the boost of the di-muon system, allow¬ 
ing separation from the signal to some extent. The processes 
e + e~ -4 vvji^fx and yy—>• v jU Vj 1J u + p~ , have slightly dif¬ 
ferent distributions of the helicity angle from the signal. All 
processes with one or two spectator electrons show signifi¬ 
cant differences from the signal, primarily in the distribution 
of the visible energy (Figure[2]i. These processes are also ef¬ 
fectively suppressed by the pj(pi) + Pt{P 2) observable. In 
addition, for the e~e —> e * e " p 1 u process, the distribu¬ 
tion of pj(pp) exhibits a peak at lower values than the sig¬ 
nal (Figure [3]). This peak corresponds to events in which the 
di-muon system recoils against electron spectators or out¬ 
going photons that are emitted below the angular cut of the 
very forward EM-shower tagging. The above is illustrated 
in Figure [3] showing the p\ distributions for representative 
background processes. 


Fig. 2 Distribution of the visible energy for the signal and e ± y—> 
e ± p + p~, e + e~ —> vvp + p~, e ± y —> e ± v^ l v fl p + p^ and e + e~ —» 
e + e~p + p~ background. 
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Fig. 3 Distribution of the pr(pp) for the signal and e + e —> 
vvp + p~, 77— > u , e~e + — > e + e~ p + p~ background with 

spectator electrons. 



P t (MM) (GeV) 


The distribution of the BDT classifier variable for the 
signal and the main background processes is shown in Fig¬ 
ure [4] (a). The classifier cut position was selected to max¬ 
imise the significance, defined as N s /y/N s +N b , where N s 
and Nb are the number of selected signal and background 
events, respectively. A plot of significance as a function of 
the position of the BDT cut is shown in Figure [4] (b). The 
optimal cut position was found at BDT = 0.23. Distributions 
of the di-muon invariant mass before and after the MVA 
selection are shown in Figure [5] Figure [5] (a) includes all 
events that pass the preselection, while Figure [5] (b) shows 
all events passing the MVA selection. All samples are nor¬ 
malised to the integrated luminosity of 1.5 ab 1 . The sig¬ 
nal preselection efficiency is 82%. The MVA selection ef¬ 
ficiency for the signal is 32%, reflecting the fact that sen¬ 
sitive observables have limited power to discriminate be¬ 
tween the signal and background. The overall signal effi¬ 
ciency including reconstruction, preselection, losses due to 
coincident tagging of Bhabha particles and the MVA is 24%, 
resulting in an expected number of 19 signal events. 


where L stands for the integrated luminosity and e s is the 
total counting efficiency for the signal, including the recon¬ 
struction, preselection and MVA selection. In the experi¬ 
ment, the number of signal events N s will be determined by 
fitting the di-muon invariant mass distribution with a func¬ 
tion f(m ^y. 

f{mpp) = N s f s (m pp ) +N b f b (m^ l ) (2) 

where f sb are probability density functions (PDF) used to 
describe the signal and the sum of all background processes, 
and N s and N b are the respective numbers of signal and back¬ 
ground events in the fitting mass window. In this analysis, 
an unbinned likelihood fit, with all parameters of f s fi{ m pp) 
fixed, is performed on simulated signal and background sam¬ 
ples. N s and N b are left as free parameters determined from 
the fit. The way the signal and background PDFs are ob¬ 
tained is discussed in Section l7TTl 

In order to estimate the statistical uncertainty of the sig¬ 
nal count, 5000 toy Monte Carlo (MC) experiments are per¬ 
formed, where pseudo-data are obtained by randomly pick¬ 
ing the signal /n uu values from the fully simulated signal 
sample, while background m pp values are randomly gen¬ 
erated from the total background PDF The size 

of the signal sample N' s and sample sizes N' b ■ of individual 
backgrounds considered, are obtained from the Poisson dis¬ 
tribution for the integrated luminosity of 1.5 ab 1 , taking 
into consideration corresponding cross-sections a and the 
selection efficiencies e {{N[) = L • o s - £ s . (N b ^ = L ■ < 7 , • £,, 
where i is indexing the different background processes listed 
in Table [TJ. 

For each toy MC experiment, the m distribution is fit¬ 
ted by the function f(m pp ) given in Eqj2] and the standard 
deviation of the resulting distribution of N s over all toy MC 
experiments is taken as the estimate of the statistical uncer¬ 
tainty of the a(Hvv) x BRt ll —> pp ) measurement. 

As will be discussed in Section [9T| the di-muon invari¬ 
ant mass distribution is sensitive to the detector pj resolu¬ 
tion, while the Higgs width I# can be considered negligible 
in comparison to the detector energy resolution. 


7 Di-muon invariant mass fit 

The quantity a(Hvv ) x BR(H —> pp ) is determined from 
the equation: 

a(Hvv) x BR(H —)• p + p~) = (1) 

L • £ s 


7.1 Signal and background PDFs 

Fully simulated samples of signal and background (Table[I]( 
are fitted to extract the PDFs. The sizes of the samples vary 
from several tens of thousands of events for the signal, up to 
a few million of events for various background processes. 

The signal PDF was defined as a linear combination of 
a Gaussian function with exponential tails, f exp and a Gaus¬ 
sian function with tails that asymptotically approach con- 
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Fig. 4 Stacked histograms of the BDT output variable for signal and background processes (a); Significance as a function of the BDT cut value 
(b). 
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Fig. 5 Stacked histograms of the di-muon invariant mass distributions with preselecton only (a) and after MVA selection (b). The distributions 
are normalised to the same integrated luminosity of 1.5 ab~*. The same legend applies to (a) and (b). 
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and 


The parameters of Eq. [3]are determined by fitting the di¬ 
muon invariant mass distribution for the signal (Figure|6]l. 

The total background PDF is defined as a linear combi¬ 
nation of a constant and exponential term: 

f b = Po-(pie P2{m - mH) + (l-pi)) (4) 

The di-muon invariant mass fit of the total background 
is shown in Figure [7] together with the fit results for the 
free parameters in Eq|4] As the normalisation to the com¬ 
mon integrated luminosity requires different normalisation 
coefficients for different processes, binned data were used 
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to combine the background processes in a straightforward 
manner and a binned fit was performed. The % 2 /Ndf of 
the background fit was 62/61. 

Fig. 6 Distribution of the di-muon invariant mass m^, of the signal 
after MVA selection, used for the PDF determination. The distribution 
is fitted with the function f s given in Eq[3] 



7.2 Distribution of the signal count 

The overall function /( m (Eqj2ji is fitted to the pseudo¬ 
data of each toy MC experiment using the unbinned likeli¬ 
hood fit. An example of a toy MC fit is given in Figure [8] 

The standard deviation of the resulting signal count dis¬ 
tribution in 5000 repeated toy MC experiments corresponds 
to the statistical uncertainty of the measurement and is 38%. 
(Figure [9]). According to EqJT] it translates into the statisti¬ 
cal uncertainty of the g(Hvv) x BR(H — » li + fl~) measure¬ 
ment, having in mind that the total uncertainty of the inte¬ 
grated luminosity can be determined at the permille level 

es. 

The statistical uncertainty of the signal counting is dom¬ 
inated by contributions from the limited signal statistics and 
from a presence of irreducible backgrounds. To estimate the 
significance of the signal against the null-hypothesis, an¬ 
other set of 5000 toy MC experiments was performed with 
zero signal count, and f(m ^) (Eq[2j) was fitted to the pseudo¬ 
data. The resulting N s distribution was centered on zero with 
a standard deviation of 5.4. Thus, in the case where the SM 
expected number of 19 signal events are found in an experi¬ 
ment, the corresponding signal significance would be 3.7 a. 

The Higgs coupling to muons, is optimally ex¬ 

tracted in a global fit procedure taking into account all Higgs 
measurements at the 350 GeV, 1.4 TeV and 3 TeV stages. 
The global fit serves to extract Higgs couplings from all 
measurements, as well as the experimental Higgs width f//. 


Fig. 7 Distribution of the di-muon invariant mass m^, of the total 
background, after MVA selection. The distribution is fitted with the 
function fb given in Eq.^] 



Fig. 8 Fitted distribution of the di-muon invariant mass for the 
sum of the signal and the total background in one toy MC experiment. 
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Because o(Hvv) x BR{H -> 8hw ^ h ^ and hav¬ 

ing access to f// and gHWW from other measurements, ex¬ 
traction of gHun is possible solely from the measurement 
presented here. An example of a minimal set of measure¬ 
ments giving a model-independent access to T// and guww 
is the following: the H -^bb measurements at both 350 GeV 
and 1.4 TeV give access to the ratio gHWW the recoil mass 
measurement at 350 GeV CM energy gives access to gnzz, 
and the H —> W 1 W measurement at 1.4 TeV gives access 

4 

to the ratio J4). The contributions of these measure¬ 

ments towards the final Agn^^ is negligible at the second 
significant digit. 











































































































Fig. 9 Distribution of the number of signal events in 5000 toy MC 
experiments. 
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The dominant contribution to the gHu/i coupling uncer¬ 
tainty is the statistical uncertainty of the measurement pre¬ 
sented here. Systematic uncertainties affect the total uncer¬ 
tainty of gHun determination only at the third significant 
digit, and thus can be neglected (Section[9]l. Under these as¬ 
sumptions, the relative uncertainty of gH/i/i is approximated 
to be 19%. 


8 Impact of electron polarization 

If 80% left-handed polarisation of the electron beam is as¬ 
sumed during the entire operation time at 1.4 Te V, the Higgs 
production cross-section through WW -fusion would be en¬ 
hanced by a factor 1.8 61- The most important background 
contribution after the MVA selection, the e + e~ —» V e V e fl p 
process, is enhanced by the same factor because it is also 
mediated by W bosons which have only left-handed interac¬ 
tions. The process e ± y—> e ± p + p~ is enhanced by a factor 
1.32, while cross-sections for other background processes 
are not significantly changed w.r.t. the unpolarised case. The 
overall selection efficiency of the signal is 30%, because the 
classifier cut position is moved to a lower value which con¬ 
sequently leads to a higher signal efficiency. The final sta¬ 
tistical uncertainty of the g(Hvv) x BR(H —► p + p~) mea¬ 
surement is 25%. The corresponding uncertainty of gH/i/i is 
13%. A summary of the results of the presented analysis is 
given in Table [3] It is important to note that all kinematic 
variables are unaffected by the beam polarization. 


9 Systematic uncertainties 

From EqJT] it is clear that uncertainties of the integrated lu¬ 
minosity and muon identification efficiency influence the 


Table 3 Summary of the cj(Hvv) x BR(H —> p + p~) measurement at 
1.4 TeV CLIC with unpolarised and 80% polarized electron beams. All 
uncertainties are statistical. 



Unpolarised 

Polarised (80%, 0%) 

N s 

19.3 ±0.1 

35 ±9 

Es 

24% 

25% 

S(a(Hvv) xBR[H^> P + p -)) 
a{Hvv)xBR{H^n+n~) 

38% 

25% 


19% 

13% 


uncertainty of the H —> p + p~ branching ratio measurement 
at the systematic level. It has been shown that at 3 TeV CLIC 
(35], where the impact of the beam-induced processes is the 
most severe, the luminosity above 75% of the nominal CM 
energy can be determined at the permille level, using low- 
angle Bhabha scattering. Below 75% of the nominal CM 
energy, the luminosity spectrum can be measured with a pre¬ 
cision of a few percent using wide-angle Bhabha scattering 
(36l . About 17% of all Higgs production events occur at a 
CM energy below 75% of the nominal CM energy. Having 
in mind the intrinsic statistical limitations of the signal sam¬ 
ple, this source of systematic uncertainty can be considered 
negligible. 

On the detector side, an important systematic effect is 
the uncertainty on the transverse momentum resolution, be¬ 
cause it directly influences the expected shape of the sig¬ 
nal m pp distribution. The sensitivity of the signal count to 
the accuracy of the knowledge of the p-\ resolution a P[ has 
been studied by performing the analysis with an artificially 
introduced uncertainty of an exaggerated magnitude on the 
assumed pj resolution used to extract the signal PDF. Re¬ 
sults of the relative shift in signal counts w.r.t. the relative 
shift of G P1 are shown in Figure[l0] The relative bias in sig¬ 
nal counting per one percent change of g pt is 0.35%. 

The uncertainty of the muon identification efficiency will 
directly influence the signal selection efficiency. In addition, 
the uncertainty of the muon polar angle resolution impacts 
the m pp reconstruction. Based on the results of the LEP ex¬ 
periments ED, it can be assumed that these detector related 
uncertainties are below a percent. 

The systematic uncertainty of the signal count caused by 
the fit with f npp defined in Eq[2j was found to be about 1% 
which is small compared to the statistical error. 

Because of the forward EM shower tagging, about 7% 
of all events are rejected by coincident detection of Bhabha 
events. This fraction must be precisely calculated taking into 
account Bhabha event distributions, beam-beam effects, as 
well as the dependence of the tagging efficiency on energy 
and angle of the incident electrons and photons. This is work 
in progress ff384]4Ql , but the uncertainty of this effect is also 
expected to be negligible compared to the statistical uncer¬ 
tainty of the measurement. 
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Fig. 10 Impact of the uncertainty of the muon pj resolution on the sig¬ 
nal counting. The relative shift of the signal count is given as a function 
of the relative shift of the pj resolution. 



Fig. 11 Dependence of the relative statistical uncertainty of the 
a(Hvv) x BR(H —» p + fi~) on the transverse momentum resolution, 
5] / px , averaged over the signal sample in the whole detector. 


cc 


X 100 - 
cc : 
m 

x 80 - 
60 y 

40 7 
20 7 
0 


CLICdp 


10 ' 


• Fast simulation 
X Full simulation 


^0 A 


Sp-i (GeV' 1 


9.1 Benefit of a improved p-\ resolution 

To estimate the benefit of a better p-\ resolution, the anal¬ 
ysis was repeated by substituting the muon four-momenta 
reconstructed in the full simulation of the signal by the four- 
momenta obtained by a parametrisation of the momentum 
resolution for several different values of the detector reso¬ 
lution. Figure o displays the approximate dependence of 
the statistical uncertainty of the measurement on the aver¬ 
age transverse momentum resolution in the whole detec¬ 
tor. Even a large improvement of the muon momentum res¬ 
olution would result in only a moderate improvement of 
the statistical uncertainty of the measured product of the 
Higgs production cross-section and the branching ratio for 
the H —> p p decay. 


10 Conclusions 

It has been shown that the measurement of the cross-section 
times the branching ratio for the SM Higgs decay into two 
muons can be performed with a relative statistical uncer¬ 
tainty of 38% at 1.4 TeV CLIC, assuming 1.5 ab 1 inte¬ 
grated luminosity with unpolarised beams. The result is dom¬ 
inated by the limited signal statistics and the irreducible back¬ 
ground. The systematic uncertainties are negligible in com¬ 
parison to the statistical one. This translates into a relative 
uncertainty of the coupling of Higgs to muons gHpp. of aprox- 
imately 19%. If the same integrated luminosity is collected 
with 80% left-handed polarisation for the electrons, the rel¬ 
ative statistical uncertainty improves to 25% and 13% for 
<j(Hvv) x BR(H —)■ p + p~) and gHpp, respectively. 
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